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BACKGROUND OF THE INVENTION 

Cross Reference to Related Applications 

This application is a continuation of application 10/167,101 filed 06/10/2002 by 
the same inventor and assigned to the same assignee. 

Field of the Invention 

The invention relates generally to GPS clocks and more particularly to an indoor 
GPS clock using long integration periods for acquiring and tracking GPS signals lower 
than about -143 dBm and providing disciplined frequency and time standard signals. 

Description of the Prior Art 

Conventional GPS positioning receivers provide GPS-based time as a byproduct 
of the resolution of a three dimensional GPS-based location. Typically, this time is 
issued in the form of a data packet that identifies the GPS-based times of the pulses of a 
one pulse per second (1 PPS) output signal. However, the accuracy for such time is 
limited by the cycle period of the local reference oscillator in the GPS receiver and any 
uncalibrated electrical length that the signal travels before it is used. There are several 
GPS timing applications where these limitations are not acceptable. Therefore, GPS 
timing applications commonly require a higher frequency signal, for example 10 MHz, 
for use as an accurate time base for maintaining an internal time standard. In order to 
meet the needs of these applications, a special type of GPS receiver, termed a GPS clock, 
has been developed. GPS clocks use supplementary techniques, such as clock bias 
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feedback, for providing an output signal having a frequency that closely tracks the 
frequency of the GPS signal. 

In some cases accurate time and/or frequency are required at a user location that is 
inside of a building. It is well-known that GPS receivers do not function well within a 
building because the building attenuates the GPS signal to a level that is too low for 
acquisition and/or tracking. It is sometimes possible to avoid this problem by placing the 
GPS signal antenna outdoors or by a window with a clear sky view and then conducting 
an amplified antenna output signal or the frequency and time signal to where it is needed. 
However, there are certain circumstances where this is impractical. 

A clear-view GPS signal available to a GPS antenna on earth from a GPS satellite 
is specified as -160 dBW (-130 dBm) by the GPS system specification ICD-GPS-200 Rev 
C published by Arinc Research Corporation of El Segundo, California published 10 
October 1993 and last revised 1 1 October 1999. Misra and Enge in "Global Positioning 
System" ISBN 0-9709544-0-9 page 288 also show -160 dBW (-130 dBm) as received 
power available to an isotropic antenna from a satellite at the zenith. Conventional GPS 
receivers acquire GPS signals as low as -130 dBm using a signal integration period 
matching the C/A PRN code epoch period of one millisecond. In addition, it is known by 
those skilled in the art that an integration period matching the GPS bit time of 20 
milliseconds can be used for achieving a processing gain of 13 dB in order to acquire 
and/or track GPS signals as low as about -143 dBm. 

The attenuation of the building will vary a great deal depending upon the type of 
building and the depth within the building. However, it is typically true that a building 
attenuates a GPS signal by more than 13 dB so that the GPS signal is lower than the 
signal acquisition that is achievable with a 20 millisecond integration. Recently, many 
workers have applied a great deal of energy to inventing techniques for receiving low 
level (that is more than 13 dB lower than open sky) signals. Unfortunately, to date these 
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techniques have not been applied in an indoor GPS clock that is capable of providing the 
frequency and time accuracies that are required. 

One of the problems that must be resolved for low GPS signal levels is the 
problem of GPS signal carrier tracking. Conventional GPS receivers employ closed loop 
carrier phase or carrier frequency feedback corrections for tracking a GPS signal. 
However, below 27 dB-Hz C/N 0 (-143 dBm referred to the antenna) and 24 dB-Hz C/N 0 
(-146 dBm referred to the antenna) typical GPS carrier phase and frequency lock loops, 
respectively, are not able to lock reliably. Existing military type GPS receivers have used 
a technique of carrier aiding by Doppler calculations. However, techniques used in 
military receivers are not directly applicable to a high performance GPS clock because a 
military GPS receiver must operate in a high dynamic environment whereas a GPS clock 
is expected to operate in an environment that is stationary. Furthermore, the GPS clock 
must provide an accurate frequency standard signal as an output that tracks the carrier 
frequency of the GPS signal. 
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SUMMARY OF THE INVENTION 



It is therefore an object of the present invention to provide an indoor GPS clock 
capable of using low level GPS signals for providing accurate GPS-based frequency and 
5 time. 

Briefly, the indoor GPS clock of the present invention is a novel combination of 
features using a processing gain of greater than 13 dB for acquiring and tracking low level 
GPS signals for providing a highly accurate GPS-based frequency. In a base coherent 

10 embodiment, the indoor GPS clock uses GPS data bit length (20 millisecond) coherent 
integrations for acquisition and tracking of low level GPS signals. In an extended 
coherent embodiment, the indoor GPS clock uses coherent integration periods longer than 
the GPS bit data for acquisition of low level GPS signals. In a coherent-incoherent 
embodiment, the indoor GPS clock incoherently combines coherent integration periods 

15 for acquisition of low level GPS signals. The low level GPS signals are then tracked with 
carrier-less tracking using GPS data bit length coherent integration periods. A clock bias 
feedback loop provides feedback for disciplining frequency and time signals. A holdover 
driver compensates for drift in the disciplined frequency and time signals for at least 
several hours in the absence of the GPS signal. 

20 

An advantage of the present invention is that an accurate GPS-based frequency is 
provided within a building where GPS signal levels are lower than about -143 dBm. 

These and other objects and features of the present invention will no doubt 
25 become obvious to those of ordinary skill in the art after having read the following 
detailed description of the preferred embodiments which are illustrated in the various 
figures. 



30 
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IN THE DRAWINGS 



Fig. 1 is a block diagram of an indoor GPS clock of the present invention; 

Figs. 2 A and 2B are block diagrams of frequency and time domain embodiments, 
respectively, of a long integration correlation machine of the indoor GPS clock of Fig. 1; 

Figs. 3 A, 3B, and 3C are block diagrams of a base coherent, an extended coherent, 
and a coherent-incoherent embodiment, respectively, of a code acquisition process of the 
indoor GPS clock of the Fig. 1; 

Fig. 4A is a timing diagram for the indoor GPS clock of Fig. 1; and 

Figs. 4B and 4C are tables for the embodiments of Figs. 3B and 3C, respectively. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 



Fig. 1 is a block diagram of an indoor GPS clock of the present invention referred 
to with a reference number 10. It is expected that the indoor GPS clock 10 is maintained 
in a single user location for any one period of operation. To initialize the operation, the 
indoor GPS clock 10 receives information for an accurate user location, recent satellite 
ephemeris parameters, and time to within ±1/2 millisecond. In an exemplary case, the 
user location can be determined with survey techniques and the ephemeris parameters can 
be received and updated from a network connection. 

The initial ±1/2 millisecond time can be determined by the indoor GPS clock 10 
itself by operating first outdoors in the manner of a conventional GPS receiver. The 
indoor GPS clock 10 is then carried indoors to its user location. Or, when one strong 
GPS signal is available, for example from a pseudolite or a GPS satellite with line of 
sight through a window, the indoor GPS clock 10 can determine the initial ±1/2 
millisecond time from that pseudolite or GPS satellite in the manner of a conventional 
GPS receiver. The initial time can also be determined outdoors with a handheld GPS 
receiver that is capable of maintaining and providing ±1/2 millisecond time for at least a 
few minutes after it loses the GPS signal when it is carried indoors and connected to the 
GPS clock 10. 

United States Patent Applications 09/860,125 entitled "Signal Receiver Using 
Coherent Integration in Interleaved Time Periods for Signal Acquisition at Low Signal 
Strength" and 09/893,179 entitled "Signal Receiver for Integrating and Combining 
Integrations in Alternating Time Segments for Signal Acquisition at Low Signal 
Strength" are incorporated herein by reference. 

The indoor GPS clock 10 includes an antenna 12, a frequency downconverter 14, 
a disciplined reference oscillator 16, a long integration correlation machine 18, a carrier 
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numerically controlled oscillator (NCO) 22, a pseudorandom (PRN) code generator 24, a 
signal processor 26, a Doppler estimator 28, a frequency and time signal synthesizer 32, 
and a navigation processor 34. The antenna 12 converts an airwave GPS signal to a 
conducted form. The frequency downconverter 14 uses a reference signal from the 
reference oscillator 16 for downconverting the radio frequency (RF) of the conducted 
GPS signal to a lower frequency and providing digital I (in-phase) and Q (quadrature 
phase) samples of the GPS signal to the correlation machine 18. 

The carrier NCO 22, the PRN generator 24, and other elements of the indoor GPS 
clock 10 use the disciplined reference signal as a time base. The correlation machine 18 
correlates the GPS signal samples with a carrier mixing signal from the carrier NCO 22 
and a PRN code replica signal from the PRN generator 24 and integrates the correlations 
for selected coherent integration periods for providing I and Q correlation data to the 
signal processor 26. 

The signal processor 26 has a signal acquisition mode and a signal tracking mode 
shown schematically with an acquisition/tracking (A/T) switch 42. Fig. 1 shows the A/T 
switch 42 in the tracking mode. The correlation machine 18 provides coherent 
integration periods for signal acquisition and/or signal tracking that are at least ten and 
preferably about twenty epochs (20 milliseconds for the C/A code) of the PRN code 
period of the GPS signal. In addition, the signal processor 26 has an extended coherent 
mode, described below, using data stripping for extending the coherent integration 
periods beyond 20 milliseconds. 

The signal processor 26 has a carrier acquisition process 44, a code acquisition 
process 45, a code tracking process 46, and a carrier hold 47. These processes 44-47 are 
implemented in a shared combination of hardware and embedded software in a digital 
signal processor (DSP), a microprocessor and a memory. The following text describes 
the processing of a single channel for a GPS signal for a single GPS satellite. However, it 
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should be understood that parallel channels may be used for processing the GPS signals 
from several GPS satellites more or less simultaneously. 

The Doppler estimator 28 uses the orbital ephemeris parameters for the satellite, 
5 time, and user position for estimating Doppler shift of the GPS signal. The carrier 

acquisition process 44 uses the estimated Doppler shift for presetting the carrier NCO 22 
to a carrier mixing frequency so that the effect of the estimated Doppler shift is 
eliminated. The code acquisition process 45 observes and processes the I (in-phase) and 
Q (quadrature phase) correlation data, described in greater detail below, for controlling 

10 the PRN code generator 24 to select a particular replica PRN code and then to sequence 
through the possible phase offsets of the code. Several PRN codes may be tried. 
Detection of the GPS signal is indicated when some function of the levels of the I and Q 
correlation data meets a certain threshold. If the GPS signal is not detected after all phase 
offsets have been tried, the carrier acquisition process 44 sets the carrier NCO 22 to 

15 another replica frequency and the code acquisition process 45 is repeated. 

When the GPS signal is detected, the carrier acquisition process 44 sets the carrier 
hold 47 so that the corresponding carrier mixing frequency is maintained by the carrier 
NCO 22. The code tracking process 46 uses the PRN code and code phase offset to 

20 initiate the signal tracking mode. In the signal tracking mode the code tracking process 
46 uses the I and Q correlation data for coherent tracking integration periods of about 
twenty epochs (20 milliseconds for GPS C/A code) to provide feedback for adjusting the 
phase offset of the PRN generator 24. This feedback is filtered in the code tracking 
process 46 in order to track GPS signals lower than -143 dBm. The Doppler estimator 28 

25 continues to compute Doppler shifts for the satellites and passes the new Doppler shifts to 
the carrier hold 47. The carrier hold 47 readjusts the carrier NCO 22 accordingly so that 
the carrier mixing frequency remains at least within an estimated ± 25 Hertz of the actual 
GPS signal. 
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It should be noted that the PRN code generator 24, the correlation machine 18 and 
the code tracking process 46 form a signal tracking loop referred to by a reference number 
48. Importantly, the indoor GPS clock 10 uses a technique described herein as carrier- 
less tracking. In conventional GPS receivers, the GPS signal is tracked with two loops, a 
code tracking loop and a carrier tracking loop. The carrier tracking loop in typical GPS 
receivers uses carrier phase or carrier frequency derived from the carrier of the GPS 
signal. In the typical GPS receiver the I and Q correlation data is used for providing a 
frequency feedback correction adjustment to a carrier NCO so that the carrier mixing 
frequency tracks the carrier frequency of the GPS signal. In the indoor GPS clock 10, 
changes in the mixing frequency from the NCO 22 are controlled by the Doppler 
estimator 28 and the carrier hold 47 and not by the I and Q correlation data. That is, the 
carrier tracking loop process that is found in typical GPS receivers and/or GPS clocks has 
been eliminated. 

The code tracking process 46 provides pseudoranges to the navigation processor 
34. The navigation processor 34 uses the pseudoranges with the user position, ephemeris 
orbital parameters, and time in a set of GPS navigation equations for resolving a clock 
bias error. The clock bias error corresponds to the phase error of the reference oscillator 
16. The navigation processor 34 also uses the pseudoranges with the user location, 
ephemeris orbital parameters, and time for providing time align information to align the 
correlation machine 18 and the signal processor 26 with the 20 millisecond transitions in 
GPS data bits. 

The indoor GPS clock 10 also includes a clock bias loop 51 and a holdover driver 
52. The clock bias loop 51 includes a clock bias filter 53 and a digital-to-analog 
converter (DAC) 54. The clock bias filter 53 receives the clock bias error from the 
navigation processor 34. The clock bias filter 53 may include differentiators, integrators, 
amplifiers and/or filters for setting the bias loop time constant, stabilizing the loop, and 
providing a conditioned clock bias error to the DAC 54. The bias loop time constant is 
between one second and four hundred seconds and preferably between one hundred fifty 
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seconds and two hundred fifty seconds. The DAC 54 converts the conditioned clock bias 
error from digital to analog form and uses the analog clock bias error to control the 
reference oscillator 16 in order to provide the disciplined reference signal. It should be 
noted that the clock bias loop 51 also includes the reference oscillator 16, the frequency 
downconverter 14, the signal tracking loop 48 (including the PRN code generator 24, the 
correlation machine 18, and the code tracking process 46), and the navigation processor 
34. 

The synthesizer 32 uses the disciplined reference signal as the time base for 
synthesizing frequency and time standard signals at various frequencies such as 10 MHz 
and 1 pulse per second (PPS). The frequency and time standard signals have a frequency 
accuracy of one part in 10" 1 1 after one day of continuous operation. The synthesizer 32 
issues the frequency and time signals as outputs from the indoor GPS clock 10. It should 
be noted that the time-of-transmission of the 1 PPS signal and the phase of the 10 MHz 
signal are synchronized with the phase of the disciplined reference signal to track GPS- 
based time and carrier phase, thereby eliminating the cycle period ambiguities of the 1 
PPS and 10 MHz signals that are found in a conventional GPS receiver. 

The indoor GPS clock 10 has a normal operation mode and a holdover mode 
shown schematically with a holdover switch 58. Fig. 1 shows the holdover switch 58 in 
the normal mode. In the normal mode, the code tracking process 46 tracks the code in the 
GPS signal; the carrier hold 47 provides updates of the Doppler shift estimates; and the 
clock bias loop filter 53 provides a closed loop correction to the disciplined reference 
oscillator 16 through the DAC 54. In the holdover mode, in the absence of GPS signal 
tracking, the holdover driver 52 provides a predicted conditioned clock bias error to the 
DAC 54 to compensate for expected temperature drift and expected aging of the reference 
oscillator 16. The aging may be initially determined when the indoor GPS clock 10 is 
manufactured or may be learned while the code tracking process 46 is tracking the GPS 
signal. The temperature drift may also be initially determined when the indoor GPS clock 
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10 is manufactured or may be learned with the aid of an internal temperature sensor while 
the code tracking process 46 is tracking the GPS signal. 

The term "coherent integration period" is used herein when GPS signal correlation 
levels are accumulated (integrated) linearly. For example, if Ii and Qi, I 2 and Q 2 , and I 3 
and Q3 are GPS signal correlations during a coherent integration period, then the coherent 
integrated (accumulated) I correlation level is Ii + I 2 + 13; the Q correlation level is Qi + 
Q2 + Q3; and the correlation power is proportional to (Ii + 1 2 + 13) 2 + (Qi + Q 2 + Q3) 2 . 
The processing gain for coherent integration increases by 3 dB for each doubling of 
length of the coherent integration period. The term "incoherent sum" is used herein when 
GPS signal correlation powers are summed. For example, if Ii and Qi, I 2 and Q 2 , and I 3 
and Q3 are GPS signal correlations during three coherent integration periods, then the 
incoherent sum of the correlation powers is proportional to Ii 2 + Qi 2 + 1 2 2 + Q 2 2 + 13 2 + 
Q3 2 . The incoherent summation increases by 1 .5 dB for each doubling of the number of 
(equal length) coherent integration periods. Therefore, coherent integration offers twice 
the processing gain. However, a factor to two longer coherent integration period tightens 
the accuracy and stability requirements of locally generated signals by a factor of two; 
whereas doubling the number of coherent integration periods in an incoherent summation 
has little or no effect on these requirements. 

Fig. 2 A is a block diagram of a time domain embodiment 18A of the long 
integration correlation machine 18. The correlation machine 18A includes I and Q carrier 
multipliers 721 and 72Q, I and Q code multipliers 741 and 74Q, and I and Q accumulators 
761 and 76Q. The I and Q carrier multipliers 721 and 72Q multiply I s and Q s GPS signal 
samples from the frequency downconverter 14 by Ic and Qc carrier mixing samples from 
the carrier NCO 22 for providing I and Q intermediate products. The I and Q code 
multipliers 741 and 74Q multiply the I and Q intermediate products by I PRN and Q PRN code 
replica samples from the PRN code generator 24 for providing I and Q code correlations. 
The I and Q accumulators 761 and 76Q use the 20 millisecond time alignment (see Fig. 
4A) from the navigation processor 34 for accumulating the I code correlations and the Q 
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code correlations, respectively, for 20 millisecond integration periods to provide the I and 
Q correlation data (I L vl and Qlvl) to the digital signal processor 26. 

Fig. 2B is a block diagram of a frequency domain embodiment 18B of the long 
integration correlation machine 18. The correlation machine 18B includes I and Q carrier 
multipliers 721 and 72Q, a signal fast Fourier transform (FFT) converter 82, a code FFT 
converter 84, 1 and Q code multipliers 861 and 86Q, a signal inverse fast Fourier 
transform (IFFT) converter 88, and I and Q accumulators 921 and 92Q. The I and Q 
carrier multipliers 721 and 72Q multiply I s and Q s GPS signal samples from the 
frequency downconverter 14 by Ic and Qc carrier mixing samples from the carrier NCO 
22 for providing I and Q intermediate products. The signal FFT converter 82 transforms 
the I and Q intermediate products into the frequency domain. The code FFT converter 84 
transforms the PRN code into the frequency domain. 

The I and Q code multipliers 861 and 86Q multiply the frequency domain I and Q 
intermediate samples by the frequency domain I PRN and Q PRN code samples to provide 
frequency domain I and Q code correlations. The signal IFFT converter 88 converts the 
frequency domain I and Q code correlations to the time domain. The I and Q 
accumulators 921 and 92Q use the 20 millisecond time alignment (see Fig. 4A) from the 
navigation processor 34 for accumulating 20 iterations of the 1 millisecond epoch of the I 
code correlations for each PRN code phase separately for providing I correlation data 
(Ilvl) and accumulating 20 iterations of the 1 millisecond epoch of the Q code 
correlations for each PRN code phase separately for providing Q correlation data (Qlvl). 

For an exemplary case of V 2 chip spacing, 1023 code phases, and 20 code epochs 
of coherent integration, the correlation machine 18 provides Ilvl correlation data for 
2*1023*20 = 40,920 phase offsets and Qlvl correlation data for 2*1023*20 = 40,920 
phase offsets. 
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There are three embodiments of the code acquisition process 45, a base coherent 
embodiment illustrated in Fig. 3 A, an extended coherent embodiment illustrated in Fig. 
3B, and a coherent-incoherent embodiment illustrated in Fig. 3C. 

Fig. 3 A is a block diagram of the base coherent embodiment 45A of the code 
acquisition process 45. The code acquisition process 45 A includes squarers 1021 and 
102Q and a summer 104. The squarers 1021 and 102Q square the I L vl and Qlvl 
correlation data. The summer 104 sums the squared levels of the I L vl and Q L vl 
correlation data for providing the sum of the squares (Ilvi 2 + Qlvl 2 )- The largest of the 
sums is compared to an acquisition threshold. When the largest sum exceeds the 
threshold, the corresponding code phase is used as the initial code phase for the code 
tracking process 46. Considering -130 dBm as a reference for the lowest GPS signal 
level that can be properly received with one millisecond coherent acquisition integration 
periods as 0 dB processing gain, the base coherent embodiment provides a processing 
gain of 13 dB, thereby enabling reception at about -143 dBm. Of course, it is understood 
by those skilled in the art that other factors may degrade performance so that the lowest 
reception levels will not be quite as good in practice as the above calculated lowest 
reception levels. It is also understood that antenna directivity can also be used 
independently to provide gain to receive lower level signals. 

Fig. 3B is a block diagram of the extended coherent embodiment 45B of the code 
acquisition process 45. The code acquisition process 45B includes an extended coherent 
processor 1 10, squarers 1 121 and 1 12Q, and a summer 1 14. The extended coherent 
processor 1 10 includes a data stripper 1 16 and accumulators 1 181 and 1 18Q. The data 
stripper 116 strips the GPS data bit information from the I L vl and Qlvl correlation data. 
For two level (one bit) information this is accomplished by inverting or non-inverting the 
Ilvl and Qlvl correlation data according to the senses of the GPS data bits for providing I 
and Q correlation data. The accumulators 1 181 and 1 18Q accumulate the stripped I L vl 
and Qlvl correlation data for providing levels of I and Q extended correlation data Ielvl 
and Qelvl for coherent integration periods of greater than the length of the data bits 
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(currently specified at 20 milliseconds in the GPS specification). For example, if the GPS 
data bits are a "1" followed by a "0", the data stripper 1 10 non-inverts the I L vl and Q LV l 
correlation data for the 20 milliseconds during the "1" data bit and inverts the I L vl and 
Qlvl correlation data for the 20 milliseconds during the "0" data bit (or vice versa). The 
senses of the GPS data bits may be obtained from a network connection. 

The squarers 1 121 and 1 12Q square the I E lvl and Qelvl correlation data. The 
summer 1 14 sums the squares of the levels of the Ielvl and Qelvl correlation data for 
providing the sum of the squares (Ielvl 2 + Qelvl 2 ). When one of the sums exceeds the 
threshold, the corresponding code phase of the largest sum is used as the initial code 
phase for the code tracking process 46. Extended coherent integration periods of 40 
milliseconds have a processing gain of 16 dB for reception as low as about -146 dBm; 
extended coherent integration periods of 100 milliseconds have a gain 20 dB for reception 
as low as about -150 dBm; and extended coherent integration periods of 400 milliseconds 
have a gain of 26 dB for reception as low as about -156 dBm. A table showing the way in 
which the extended coherent code acquisition process 45B processes the I L vl and Qlvl 
correlation data is illustrated in Fig. 4B and described in the accompanying detailed 
description below. 

For a 20 millisecond coherent acquisition integration periods the carrier mixing 
frequencies are stepped in about 25 Hertz steps; for a 40 millisecond extended coherent 
acquisition integration periods the carrier mixing frequencies are stepped in about 12 1/2 
Hertz steps; and so on where the step size is calculated in Hertz as 25 times 20 divided by 
the coherent integration period in milliseconds. It should be noted that the coherent 
acquisition integration periods are not required to be integer multiples of 20 milliseconds. 
For example, for 50 millisecond coherent integration periods, the processing gain is step 
size is preferably about 10 Hertz. 

Fig. 3C is a block diagram of the coherent-incoherent embodiment 45C of the 
code acquisition process 45 referred to by a reference number 45C. The code acquisition 

A1298CON1 "Indoor GPS Clock" by Yu 

14 



process 45C includes squarers 1221 and 122Q, incoherent summers 1241 and 124Q, and a 
summer 126. The squarers 1221 and 122Q square the levels of the I L vl and Q LV l 
correlation data for providing I LVL 2 and Qlvl 2 for "N" of the 20 millisecond coherent 
acquisition integration periods for each code phase. The incoherent summers 1241 and 
124Q sum the squares Ilvl 2 and Qlvl 2 of the correlation data Ilvl and Qlvl for each 
coherent acquisition integration period for each code phase. The summer 126 sums 
squared correlation data for providing an incoherent sum. A table showing the way in 
which the coherent incoherent code acquisition process 45C uses the I L vl and Qlvl 
correlation data is illustrated in Fig. 4C and described in the accompanying detailed 
description below. 

A coherent-incoherent period of one second made up of fifty incoherent sums of 
the 20 millisecond coherent integrations gives a processing gain of about 8.5 dB for the 
fifty incoherent periods plus 13 dB for the 20 millisecond coherent integration period for 
a total processing gain of about 21.5 dB resulting in an improvement of initial acquisition 
(no GPS signal is being tracked before acquisition) power level from -130 dBm to about 
-151.5 dBm. A coherent-incoherent period of 80 milliseconds made up of four 
incoherent sums of the 20 millisecond coherent integrations gives a processing gain of 
about 3 dB for the twenty incoherent periods plus 13 dB for the 20 millisecond coherent 
integration period for a total processing gain of about 16 dB resulting in an improvement 
of initial acquisition power level from -130 dBm to about -146 dBm. A coherent- 
incoherent period of 10 seconds made up of five hundred incoherent sums of the 20 
millisecond coherent integrations gives a processing gain of about 13.5 dB for the five 
hundred incoherent periods plus 13 dB for the 20 millisecond coherent integration period 
for a total processing gain of 26.5 dB resulting in an improvement of initial acquisition 
power level from -130 dBm to about -156.5 dBm. 

Processing gain and improvement in acquisition and tracking signal levels can be 
tested by inserting a variable attenuator in a cable path between a passive patch antenna 
and a receiver preamplifier. The receiver preamplifier may be a part of the frequency 
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downconverter 14. A reference initial acquisition and/or tracking attenuator setting is 
found for a conventional GPS receiver using one millisecond coherent integration periods 
for the lowest signal that can be acquired or tracked. The reference attenuator setting 
establishes a reference 0 dB processing gain. The improvement of processing gain with 
the present invention is then found by determining the maximum additional attenuation 
that can be inserted for acquisition and/or tracking. A GPS satellite simulator such as a 
NT Northern Telecom STR2760 GPS Satellite Simulator can also be used. 

Fig. 4A is a timing diagram of the operation of the long integration correlation 
machine 18. The diagram shows a data bit stream of GPS data bits and 1st, 2nd, 3rd, and 
so on to an Nth 20 millisecond integration period. The data bits are known to be 20 
milliseconds long. As described above the correlation machine 18 uses the time 
alignment information from the navigation processor 34 for aligning its integration 
periods to the data bits initially within ±1/2 millisecond and then more accurately after 
the GPS signal has been tracked. 

Fig. 4B is a table illustrating the way in which the extended coherent code 
acquisition process 45B processes the I L vl and Qlvl correlation data. The code 
acquisition process 45B receives I L vl and Q LV l correlation data for coherent acquisition 
integration periods (shown with subscripts 1, 2 and so on to N) and K possible code 
phases (shown with subscripts <|>1, <|>2, and so on to K). The I correlation data is Ii^i, l 2 $u 
and so on to I N ,<|>i for code phase one I\^ 2i l 2 $ 2 , and so on to I N ,<|> 2 for code phase two 
((J) 2 ); continuing to Ii^k, and so on to In^k for code phase K ((|)k). The I E lvl 
extended correlation level is (Ii^i + 1 2 ^\ + . . . + In^i) for code phase one (I\^ 2 + h#2 
+ • • • + h^i) for code phase two (fa); and so on to + I 2 ^ K + . . . + In,4>k) for code 
phase K (<t> K ). Similarly, the Q correlation data is Qi^i, Q 2 $u and so on to Qn,*i; Qi,<p, 

and so on to Q N>+2 ; continuing to Q lf+K , Q 2 ,*k, and so on to Q Ni+K . The Qelvl 
extended correlation level is (Q 1>+1 + Q 2)(j) , + . . . + Q N)<M ); (Q142 + Q2,*2 + - . - + CW); 
and so on to (Qi j4>k + Qi^k + . . . + Qn.^k). 
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The code acquisition process 45B determines the sum of the squares of the sums 
of the I and Q correlation levels [(h M + I 2 ^i + . . . + I N ,*i) 2 + CQu+i + Q241 + . . . + 
Qn,4>i) 2 ] for code phiase (<|>i); the sum of the squares of the sums [(I\$2 + h$2 + . . . + W2) 2 
+ (Q 1 ,<t>2 + Q242 + . . . + Qn^) 2 ] for code phase (<|) 2 ); and so on to the sum of the squares of 
the sums [(I UK + h#L + ■ ■ - + Wk) 2 + (Qi,*k + Q24K + - . - + Qn,*k) 2 ] for code phase (<|>k). 
The largest of the sums of squares of the sums is compared to an acquisition threshold. 
When the largest exceeds the threshold, the corresponding code phase is used as the 
initial code phase for the code tracking process 46. 

Fig. 4C is a table illustrating the way in which the coherent-incoherent code 

acquisition process 45C processes the I L vl and Q LV l correlation data. The code 

acquisition process 45C receives I L vl and Q L vl correlation data for coherent acquisition 

integration periods (shown with subscripts 1, 2 and so on to N) and K possible code 

phases (shown with subscripts <|>1, $2, and so on to K). The I correlation data is 1 2 ,<t>i, 

and so on to I N ,<|>i for code phase one I lj+2 , 1242, and so on to I N ,^2 for code phase two 

(<|> 2 ); continuing to I^k, I24K, and so on to I NiWC for code phase K (<|> K ). Similarly, the Q 

correlation data is Qi^i, Q241, and so on to Qn,*i; Qi,<f>2, Q242, and so on to Qn,^; 

continuing to Q\$k, Q24K, and so on to Qn^k. The code acquisition process 45C 

determines the sum of the squares of the I and Q correlation levels (l 2 \^\ + + . . . + 

I 2 N,<j>i + Q 2 i,<t>i + Q 2 2,<m + . . . + Q 2 n,(|>i) for code phase (M; the sum of the squares (l 2 \$ 2 + 
2222 2 

I 2,<t>2 + . . . + 1 N,*2 + Q 1,42 + Q 2,<t>2 + . . . + Q N,<p) for code phase (<|) 2 ); and so on to the 
sum of the squares (I 2 ^ + l\$ K + . . . + l\$ K + Q 2 UIC + Q 2 2 ^ K + . . . + Q 2 n,<|>k) for code 
phase (<|>ic). The largest of the sum of the squares is compared to an acquisition threshold. 
When the largest sum exceeds the threshold, the corresponding code phase is used as the 
initial code phase for the code tracking process 46. 

It will be appreciated by those having ordinary knowledge of the art that the code 
acquisition process 45 also includes search algorithms, filters and other conditioning 
elements for processing the I and Q correlation data in order to control the PRN code 
generator 24 for GPS signal acquisition. 
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Returning to Fig. 1, the correlation machine 18 provides sequential coherent 
integrations for 20 millisecond periods for prompt I correlation Ip, early minus late I 
correlation Ie-l, prompt Q correlation Q P , and early minus late Q correlation Qe-l- The 
code tracking process 46 computes 9 E rr according to an equation 1 . 

6 E rr = (l P *l E . L +Qp*Q E - L )/(ip 2 +Qp 2 ) (1) 

The code tracking process 46 filters the 0err with a narrow bandwidth, preferably about 
one Hertz, for the code adjustment feedback to the PRN code generator 24 in order to 
continue to track low level GPS signals after they have been acquired. Those having 
ordinary knowledge in the art will understand that the code tracking process 46 also 
includes filtering, carrier aiding, and other conditioning elements for minimizing noise 
and stabilizing the signal tracking loop 48. 

The indoor GPS clock 10 as described above can be implemented with a 
combination of DSP hardware and software and microprocessor hardware and software 
where the software is stored in memory devices of various types. Those skilled in the art 
will recognize the separate functions of the elements as described above can operate with 
shared hardware and be allocated in various ways between DSP devices, microprocessor 
devices, and memory. 

Although the present invention has been described in terms of the presently 
preferred embodiments, it is to be understood that such disclosure is not to be interpreted 
as limiting. Various alterations and modifications will no doubt become apparent to 
those skilled in the art after having read the above disclosure. Accordingly, it is intended 
that the appended claims be interpreted as covering all alterations and modifications as 
fall within the true spirit and scope of the invention. 

What is claimed is: 
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